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SIGNIFICANCE  AND  EXPLANATION 


Studies  employing  both  asymptotic  and  numerical  methods  have  shown  that, 
in  the  flow  of  two  fluids  with  different  viscosities,  interfacial 
instabilities  can  arise  at  any  Reynolds  number.  Among  the  geometries  that 
have  been  studied  are  plane  Couette  and  Foiseuille  flow,  pipe  flow  and  flow 
between  rotating  cylinders. 

Hie  question  thus  arises  what  happens  in  those  cases  where  the  interface 
is  unstable.  Which  new  flow  patterns  that  would  be  stable  can  develop?  C.  S. 
Yih  has  conjectured  that  there  may  be  traveling  interfacial  waves.  More 
recent  analysis  by  Hooper  and  Boyd  has  shown  that  there  is  a  crucial 
difference  between  the  cases  of  zero  and  non-zero  surface  tension  between  the 
two  fluids.  If  the  surface  tension  is  zero,  then  waves  of  arbitrarily  short 
wave  length  are  unstable.  In  this  case,  standard  methods  of  bifurcation 
theory  are  not  applicable,  and  it  seems  possible  that  no  smooth  interface, 
steady  or  unsteady,  would  be  stable.  In  fact,  we  believe  that  this  might  be  a 
mechanism  for  the  formation  of  emulsions. 

If,  on  the  other  hand,  surface  tension  is  non-zero,  then  bifurcation 
theory  is  applicable  and  Yih’s  conjecture  could  be  true  if  the  bifurcation 
turns  out  to  be  supercritical.  In  this  paper,  we  study  plane  Couette  flow 
with  periodic  boundary  conditions  imposed  in  the  streamwise  direction.  We 
consider  a  "basic"  solution  with  a  flat  interface  parallel  to  the  walls.  We 
prove  the  existence  of  bifurcating  periodic  solutions  (interfacial  waves)  at  a 
point  where  the  flat  interface  becomes  unstable. 
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HOPF  BIFURCATION  IN  TWO-COMPONENT  FLOW 


M.  Renardy^-  and  D.  D.  Joseph2 

1.  Introduction 

The  stability  of  two-component  parallel  shear  flows  has  been  analyzed  by  long-wave 
asymptotics  [5],  [16],  short-wave  asymptotics  [6],  [13]  and  numerically  [11],  (13).  These 
studies  show  that.  If  the  fluids  have  different  viscosities,  then  instabilities  can  arise 
at  all  Reynolds  numbers. 

This  raises  the  question  of  possible  alternative  flow  patterns  which  might  be 
stable.  Ylh  [16]  has  conjectured  that  wavy  interfaces  might  develop.  The  analysis  of 
Hooper  and  Boyd  [6]  reveals  a  crucial  difference  between  the  cases  of  zero  and  non-zero 
surface  tension  between  the  fluids.  If  the  surface  tension  Is  zero,  then  sufficiently 
short  waves  are  always  unstable,  l.e.,  there  Is  an  Infinite  number  of  unstable  modes. 

This  situation  Is  very  much  unlike  the  usual  problems  of  bifurcation  theory,  and  we 
believe  It  is  possible  that  no  smooth  Interface,  steady  or  unsteady,  would  be  stable  in 
this  situation.  (In  reality,  of  course,  the  surface  tension  Is  not  zero,  but  there  will 
be  instability  for  very  short  waves  when  Che  surface  tension  Is  small  and  Che  Reynolds 
number  Is  large.  We  think  that  this  Instability  mechanism  may  be  relevant  In  the 
formation  of  emulsions.) 

In  the  case  of  non-zero  surface  tension,  however,  one  can  establish  a  bifurcation 
theorem.  If  the  bifurcation  turns  out  supercritical,  this  provides  a  basis  for  Ylh's 
conjecture.  For  the  sake  of  simplicity,  we  confine  attention  to  plane  Couette  flow,  but 
It  Is  clear  that  similar  techniques  can  be  applied  to  more  complicated  geometries  such  as 
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concentric  flow  in  pipes  or  between  rotating  cylinders.  We  consider  plane  Couette  flow  of 
two  fluids  with  equal  density,  but  different  viscosities,  and  an  interface  parallel  to  the 
plates.  Periodic  boundary  conditions  are  inposed  in  the  strearawise  direction.  Evidently, 
this  configuration  is  stable  at  rest.  If  there  is  a  flow,  however,  then  instabilities  can 
develop  [6],  [16].  Since,  however,  surface  tension  will  damp  short  waves,  there  can 
always  be  only  a  finite  number  of  unstable  modes.  Generically,  as  the  flow  rate  is 
increased,  one  specific  node  will  be  first  to  become  unstable.  Since  the  eigenvalues  are 
ccaplex,  one  expects  a  bifurcation  of  the  Hopf  type  [7],  leading  to  traveling  interfacial 
waves . 

the  Hopf  bifurcation  theorem  in  infinite  dimensions  [4],  [8-10],  [14]  relies  on 
coercive  estimates  for  the  linearized  equations.  For  one-component  free  surface  flows 
such  estimates  were  derived  by  Beale  [2],  [3],  and  we  shall,  in  Chapter  3,  derive 
analogous  estimates  for  two-component  flow.  Our  proof  differs  from  Beale's  and  is 
slightly  simpler.  Using  these  coercive  estimates,  we  can  then  establish  a  bifurcation 
theorem  in  Chapter  4.  In  Chapters  5-8,  we  outline  an  algorithm  for  the  computation  of 
bifurcating  solutions. 

2.  Formulation  of  the  Problem 

We  consider  two-dimensional  flow  of  two  fluids  with  different  viscosities  and  equal 
densities  between  parallel  plates. 


♦  V 


The  motion  in  each  fluid  is  described  by  the  Navier-Stokes  equations: 

(2.1)  P(u  ♦  (u  •  V)  u)  -  h^u  -  ’p,  ^ 

V-u  -  0  T  j  0  <  y  <  h(x). 


-2- 


(2.2)  P(v  ♦  (v«7)v)  »  n2  4v  -  7q, 


h(x)  <  y  <  1. 


7*v  -  0  ,  J 

We  have  no  slip  conditions  at  the  walls: 

(2.3)  u  -  0  at  y  -  0, 

(2.4) 

v  -  (V  0)  at  y  -  1. 

—  Of 

Across  the  interface,  there  must  be  continuity  of  velocity  and  traction.  This  leads  to 
the  equations 

(2.5)  ji  ■  at  y  -  h(x), 

”  '2  JU2  3ul  •  3u2  3ul 

(2.6)  (d-h  )  (JJ-  ♦  )  +  2h  <3—  -  57-  )  > 


-2.  .  3v2  3v1 


,  3v2  3v, 


n2  {<,"h  ’  (  IST  +  ST  ’  +  ^  (  ST  '  ST  >}  at  y“  h(x)' 


3u2 

1  3y“  " 

.  3u1 
p-2hV?ir*Jsr>+h 

3v2  _  .  3vJ  . 

J— -q-2hni,^7-  +  5—  )th  (2n2  -  q) 


(iW/2 


at  y  -  h(x) . 


Here  T  is  the  surface  tension  parameter.  Finally,  we  have  the  kinematic  boundary 
condition 

(2.8)  h  +  u1  h  -  u2_ 

He  are  interested  in  solutions  to  (2.1)  -  (2.8)  which  have  a  given  period  L  in  the  x- 
direction  and  are  periodic  in  t.  We  denote  by  (1^  the  set  {(x,y)j  0  <  x  <  L, 

0  <  y  <  h(x)},  by  the  set  {(x,y)|  0  <  x  <  t,  h(x)  <  y  «  1}  and  by 

r  the  Interface  {(x,y)|  0  <  x  <  L,  y  -  h(x)}.  The  spaces  hNs^),  Hk(fl2),  Hk(T) 

consist  of  those  functions  which  have  k  square  integrable  derivatives  and  satisfy 
periodic  boundary  conditions  in  the  x-direction. 
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3.  The  linearized  problem 

In  this  chapter,  we  obtain  coercive  estimates  on  the  linear  problem,  which  we  shall 


need  later.  He  put  VQ  -  0  and  linearise  (2.1)  -  (2.8)  at  the  rest  state  =  0, 


with  a  flat  interface  h  «  hQ.  We  include  inhomogeneous  terms  in  (2.1),  (2.2),  (2.6)  and 


(2.7).  is  replaced  by  X.  This  leads  to  the  problem: 

(3.1)  Xpu  »  n  Au-^p+f  ^ 


1,  A«  “  7P  +  I,» 

f*tt  »  0  ,  J 


0  <  y  <  h 


(3.2) 

Xpv  »  hjAv  -  Vq  +  f 
V»v  -  0  , 

}  ^ 

(3.3) 

u.  “  0  ,  y  -  0, 

(3.4) 

v.  -  0  ,  y  -  1, 

(3.5) 

u  -  v  ,  y  -  hQ, 

3u  3u 

3v  3v 

(3.6) 

ni(  37"  +  3^"  *  "  n2  ( 

35T+37“ 

3u  3v 

(3.7) 

2T11  3~  “  P  “  2n2  3~ 

+  q  -  Th  * 

(3.8) 

Xh  ■  u  ,  y  m  h  . 

2  o 

He  seek  solutions  periodic  in  x  with  period  L,  and  of  course  the  same  periodicity  is 
assumed  for  f tj  and  f^.  Our  goal  is  the  following  estimate: 


Theorem  3. 1 i 

Let  y  >  0.  Por  He  X  >  Y,  the  following  estimate  holds  for  solutions  of  (3.1)  -  (3.8). 


I  u  I  -  +  I  v  1  +  Ipl  +  Iql  | X |  I  u  I 

«  <n4>  m  <n2>  H’(ot)  h'(b2>  “  lz 


(3.9)  +  M  *  1  '  ,  ♦  «hl  ,,,  +  M  Ihl  <  C  [if  I  ♦  If  I 

L  <«2)  H5/2(D  H3/Z(D  LZ(£l1)  2  L2(S)2) 


♦  'V  1/2  *  IM1/4'*3«  2  ♦  1/2  1 

J  H v  (D  3  L  (D  HVZ(D 


>  A  .V  a 
• .  ■ • . ' 


We  multiply  (3.1)  by  u  (the  complex  conjugate  of  ji),  and  (3.2)  by  v,  add  them  and 


integrate  over  the  domain.  Integrating  by  parts,  and  using  the  boundary  and  interface 


conditions,  we  find 


Xp  ( lul  2  +  ivi  2)  +1/2n1  <  u,  u  >  +  V2n2  <  v,  v  >  +  T  T  Ih 


L  L 


(3.10) 


L  L 

(£ ,.  U)  2  +  (f2  ,  v)  2  ♦  I  f3u1(y“h0>  ^  +  /  f4  u2  (y  "  hoJ  ^ 
L  L  O  o 


Here  <  u,  u  >  «  J(Vu  ♦  (Vu)T)  :  (Vu  +  (Vu)T). 


Since  u^  is  divergence  free,  u2  has  a  trace  on  the  interface  whose  H  '2  -  norm  can  be 


bounded  by  lul  [15],  hence  the  last  term  on  the  right  side  can  be  bounded  by 

“  LZ(J11  > 

If  I  lul  .  The  third  term  is  bounded  by 

4  HV2(D  ~  X.2 ( £1 , ) 


lf,l  ,  '  u  •  ,  <  ,  '  u  •  \  •  U  |  2 

3  tV)  L  (T)  J  lV)  iT^)  H^O,) 


We  assume  that  the  right  side  of  (3.9)  is  bounded  by  a  constant  of  order  one,  and  we  wish 


to  bound  the  left  side.  Prom  (3.10),  we  immediately  obtain  bounds  for 


lul  i'  ’i*  ,hl  i  < we  have  used  Korn's  inequality  here). 

H  ~  H  H 

In  the  following,  we  make  repeated  use  of  the  following  estimates: 


(3.11) 


lul  ,  +  'vl  +  Ipl  ♦  Iql 

H  H  H  H 


<  Cdf,.  t  If.  t  lf3.  ♦  |X|.h.  ♦  IM.U.  *  |X|.V.  ), 

L  la  H  H  la  la 


(3.12)  Ihl  ...  <  C(lul  ,  +  Ivl  ,  +  Ipl  ,  ♦  Iql  ,  *  If. I  _._), 
h5/2  H2  H2  H1  H1  4  HV2 


.b.  <  ih.  f3  .hi ;;; . 

H  n  H 


*.  •*.  i-. 


•»:  •  v 


/ ./  ./  .t*  •  . 


.  s'*  s'  .*• 
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->n/in 

(3.19)  e<C  (1  +  i>1  *  1 X j 


2n+1/3n+1+I  ,*^.,1/3,. 


From  this,  (3.17)  follows  easily. 


Next,  we  wish  to  show  that  (3.14)n+i  -  (3.  16)n+1  follow  from  (3.17)n.  Using  (3.8),  the 


trace  theorem  and  the  convexity  property  of  Sobolev  norms,  we  find 
C 


(3.20)  Ihl  ^  <  -r^y  /lul  ^  /lul 
H  H  H 

Moreover  (3.10)  implies 


(3.21) 

Ivl  , 

+  lul 

<  C(lul 

+  Ivl 

'  H1 

"h1 

"l2 

: 

and 

(3.22) 

'-V2 

♦  lul 

L 

<  C  [ Ihl 

H1 

*  i»i 

With 

X  = 

lul  + 

~  H 

Ivl  ,  y  = 
H 

A2 

1-1/4 


1/2,  1/2 


■1/2, 


, -3/4  1/2  1/2  ,1/2 

I  '“I  2  '-%•  ’ 


(3.22)  take  the  following  form 


(3.23)  x<  C(yV2+  |X|"V8  x  1/4  y  V4>, 


2n_1/3n-  ~ 

(3.24)  y«C(|X|  4x  1/2.  lll-V2„  1/2 


+  |X|  -  y  + 

1  /2 

From  (3.23)  it  follows  that  x  <  Cy  or  X  <  C | X  |  y 

first  or  second  term  on  the  right  is  bigger.  In  the  first  case, 
(3.24)  yields 


5/8  1/4  1/4, 

y  x  )  . 


depending  on  whether  the 


(3.25)  y<0  (IXI2"  1/3n_  V/4+  |X|-1/2y  V* 


♦  IM'5/8  y  3/8> 


From  this  (3.15)n+.j  is  immediate.  In  the  second  case,  we  get 


ll.2<|  y  <  ‘  ‘  ’2  y1'6  .  |M-’/2  y  •  IM 


-J/<-  2.  y  VI, 


From  this,  (3.15>n+1  is  also  immediate.  From  (3.15)n+1  and  (3.21)  follows  (3.16)n+1,  and 
using  (3.20)  and  <3.17)n  we  find  (3.14)n+,. 
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a 


A.N.'.A  \  -  *  *  •  «  •  •  -  -  •  .  A  v\‘  w*  -*  V  “  ‘  -  •  •  ,  •  .  ■ 


.■  ,v; 

*  j  *.»  '.ar  k  ’  <  *  •  »  »  i  w  *  . 


N  .‘.-.V.-vl 

V.v-.VVV 


A  A 


To  proceed  further,  we  take  difference  quotients  in  the  x  -  direction.  These  satisfy  the 


sane  equations  (3.1)  -  (3.8)  with  the  f's  replaced  by  their  difference  quotients.  From 

1 


(3.10),  we  then  see  that  the  H  -norms  of  all  x-der ivatives  of  u  and  v.  can  he  estimated  by 
terns  of  order  1.  The  divergence  condition  now  yields 


•u2'2,  lv2l2  <  C. 


.  i  -5/4+e 

Equation  (3.20)  now  yields  Ihl  <  O’  |  X  |  ,  and  (3.21),  (3.22)  now  imply 

C  H 

'h'2  +  's',  4  TTT  .  The  rest  follows  easily.  It  remains  to  prove  (3.11).  For  this, 
L  L 


we  only  need  to  show  that  (3.1),  (3.2)  with  boundary  conditions  (3.3)  -  (3.6)  and  (3.8) 
frost  an  elliptic  system  in  the  sense  of  Agmon,  Douglis  and  Nirenberg  [1J.  We  can  formally 
regard  (3.1),  (3.2)  as  being  posed  in  the  same  domain  by  mapping  the  strip  occupied  by 
fluid  2  onto  the  strip  occupied  by  fluid  1.  We  do  this  in  such  a  way  that  the  interface 
is  upped  onto  itself  and  the  solid  boundaries  are  mapped  onto  each  other.  This  yields  a 
system  for  the  six  unknowns  ( u ^ ,  p,  v ^ ,  v2,  q),  which  are  now  defined  on  the  same 
domain.  That  the  Stokes  equations  are  an  elliptic  system  is  well  known.  The  same  holds 
of  course  for  two  sets  of  Stokes  equations.  It  is  also  well  known  that  Dirichlet  boundary 
conditions  satisfy  the  complementing  condition.  Showing  the  complementing  nature  of  the 
interface  condition  is  a  straight  forward  calculation,  which  we  omit. 


Remarks: 

We  have  so  far  only  given  estimates  for  a  solution  that  was  assumed  to  exist  and  have 
the  regularity  implied  by  the  left  hand  side  of  (3.9).  Such  estimates  show  that  (3.1)  - 
(3.8)  for  Re  X  >  T  >  0  is  solvable  for  a  closed  set  of  f's  (in  the  topology  indicated  by 
the  right  side  of  (3.9)).  Solvability  for  a  dense  set  of  f's  can  be  shown  in  a  straight 
forward  manner  by  separation  of  variables.  From  this  we  see  that  in  fact,  for  any 
X  with  Re  X  >  0,  (3.1)  -  (3.8)  has  a  unique  solution.  (3.9)  holds  uniformly  in  any 

closed  subset  of  any  right  half-plane,  if  this  subset  contains  no  eigenvalues.  Moreover, 
by  compactness,  the  number  of  eigenvalues  is  countable,  and  there  can  only  be  finitely 
many  in  any  bounded  set.  (3.10)  implies  that  all  eigenvalues  have  negative  real  parts. 


4.  Bifurcation  to  travelling  waves 


It  is  convenient  to  use  a  domain  mapping  which  takes  the  domain  occupied  by  each  fluid  to 


a  fixed  one.  The  most  straight  forward  way  to  construct  such  a  mapping  is  to  stretch  or 


contract  vertical  lines.  We  shall  in  addition  transform  the  velocity  fields  in  such  a  way 


that  the  divergence  condition  is  preserved  and  (2.8)  reduces  the  linearized  form  even  in 


the  nonlinear  case.  In  doing  this,  we  essentially  follow  Beale  (2). 


Let  y  «  hQ  be  the  flat  interface  of  the  rest  state,  and  let  y  =  h(x,t)  be  the  actual 


interface,  which  we  assume  periodic  in  x  with  period  L. 


Let  P  be  any  linear  extension  operator  that  maps  functions  h(?)  e  H  ,  into  functions 


h(C,n)  such  that  h(C,h  )  -  h(5)  and  h  e  H  2  (there  are  many  ways  to  construct  P). 

o 


For  simplicity,  we  also  assume  that  P  takes  h  5  hQ  to  hi  hQ.  Let  then 


f  (n)  be  a  C  -  function  of  n  such  that  f  »  1  near  n  «  h  and  f„  «  0  near 
o  o  o  ° 


n  ■  0  and  n  ■  1.  Define  h(C,n,t)  »  h(C,n,t)  •  f  (n)  +  h  ( 1— f  ( n ) ) 

o  o  o 


We  now  define 


G(c,n,t)  -  (c,n  •  h(y <""> 
o 


Evidently,  0  maps  the  strip  0  <  n  <  h  to  fl  ,  and  the  strip  h  <  n  <  1  to  n  . 

O  I  O  2 


The  velocities  are  transformed  as  follows: 

30 .  . 


(4.2)  u  <fl(C,n,t)) 


3^  7 


where  J  is  the  jacobian  of  0.  Of  course,  in  SI  is  transformed  in  the  same  way. 


Explicitly,  (4.2)  reads 


(4.3)  u  «  (—  +  n  —  )  (u  ,  n  u  +  (r—  +  n  —  )  u  ). 

—  h  h  1  h  1  h  h  2 

o  o  o  o  o 


Formula  (4.2)  is  set  up  in  such  a  way  that  u,  1.  are  divergence-free  in  the  x,y-plane,  if 


u,  v  are  divergence-free  in  the  C,  n-plane.  Moreover,  (2.8)  now  assumes  the  simple 
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(4.4)  h  <5,t)  =  u  (?,h  ,t). 

z  o 

It  is  also  clear  that  (2.5)  does  not  change  its  form,  i.e.  at  y  «  h(x)  c 

becomes  u  m  v  at  n  *  h  .  The  boundary  conditions  at  the  walls  are  also  presei 
—  —  o 

When  these  substitutions  are  inserted  into  (2.1)  -  (2.8),  we  obtain  a  new  set  of 
equations,  which  we  do  not  write  down  explicitly.  We  shall  refer  to  these  new  eqvj 
as  (2.1)*  -  (2.8)*.  If  we  have  a  flat  interface  h  =  hQ,  then  of  course  (4.3)  reads 
u  »  u,  and  (2.1)*  -  (2.8)*  have  the  same  form  as  (2.1)  -  (2.8). 

Plane  Couette  flow  is  the  following  solution  of  (2.1)  -  (2.8): 


*  _ ,  *  Vo  *  *  Vo  +  ww 

=  Ul  “  VVVV  y'  “2  =  V1  =  v^Tv"  ' 


v2  ■  0,  p  =  0,  q  -  0. 

We  can  linearize  at  this  solution,  and  obtain  a  set  of  linearized  equations  analog 
(3.1)  -  (3.8).  As  usual,  we  call  X  an  eigenvalue,  if  the  homogeneous  linearized 
has  nontrivial  solutions.  The  estimates  in  chapter  3  imply  that,  for  VQ  “  0(rest 
is  a  countable  sequence  of  eigenvalues,  all  in  a  sector  of  the  left  half  plane  bou 
away  from  the  imaginary  axis.  All  these  eigenvalues  have  finite  multiplicity  and 
index  zero.  Estimates  like  those  in  Chapter  3  can  easily  be  extended  to  the  linea 
at  Couette  flow  with  finite  V0.  If  we  linearize  (2.1)  -  (2.8)  at  this  flow,  then 
are  a  number  of  terms  perturbing  (3.1)  -  (3.8).  All  these  terms  are  relatively  cc 
except  the  one  resulting  from  Uj  h  in  (2.8).  This  latter  term  vanishes  in  a  ft 
moving  with  the  fluid  on  the  interface.  Standard  perturbation  theory  [12]  can  now 
to  show  that  estimates  like  in  Chapter  3  hold  for  X  in  any  closed  set  that  lies 
right  half  plane  and  contains  no  eigenvalues.  However,  there  can,  and  as  [6],  [1< 
there  will  be  a  finite  number  of  eigenvalues  with  positive  real  parts  if  Vc  is  ] 
enough.  Generically,  there  will  be  a  critical  value  V^,  such  that,  for  Vo  <  V( 
eigenvalues  have  negative  real  parts,  but  a  pair  of  simple  complex  conjugate 
eigenvalues  crosses  the  imaginary  axis  transversally  as  V  increases  past  V  . 


denote  these  imaginary  eigenvalues  by  ±iu0* 

We  introduce  the  following  substitutions  in  (2.1) 


(2.8)*: 


T  -  «t,  u*  “  u  -  u(V  ),  v*  ”  v  -  v(V  ),  h*  ”  (h-h  )  — 

—  —  —  o  —  —  —  O  OU) 

k° 

We  use  the  following  notation  for  function  space:  H  denotes  the  space  of  functions 

defined  on  -*•<£<«,  0<n<hQ,  which  have  period  L  in  C  and  H^-regularity . 

k  k 

Similarly,  we  define  H  .  Finally  H.  is  the  set  of  k  times  differentiable  periodic 

n  r 

functions  depending  on  "  alone. 

H^CX)  denotes  the  spaces  of  all  2v-periodic  functions  defined  on  -“  <  x  <  “ ,  taking 

values  in  X,  and  having  k  square  inteqrable  derivatives. 

For  the  analysis  of  (2.1)*  -  (2.6)*,  we  choose  the  following  space  V: 

v  -  { ( u* , v* ,  P,  q,  h* ) |  U*  e  h\h~  )  nH2(ii  ),  V*  e  h’(h~  )  n  h2<iJ  ), 

1  1  2  2 

p  e  h\h~  ),  q  e  ),  h*  e  h’(h5^2)  n  h2(h3^2);  div  u*  -  o,  div  v*  -  o. 


//  p  +  //  q  «  0,  u*  -  0  at  n  -  0,  v*  -  0  at  n  “  1,  u*  “  v *  at 
“l  «2 

n  -  hQ,  u)o  -  u*  at  n  -  hQ,  /  h*(C)dC  -  0} 


Functions  in  this  space  have  sufficient  regularity  such  that  all  the  nonlinearities  in 
(2.1)*  -  (2.8)*  are  defined.  TSie  coercive  estimates  of  Chapter  3  and  the  implicit 
function  theorem  can  be  used  to  prove  the  following  result  (the  details  of  the  proof  are 
analogous  to  the  standard  proof  of  the  Hopf  bifurcation  theorem  and  are  omitted  here). 


Theorem  4. 1 : 

Assume  that,  at  VQ  -  V00,  there  is  a  pair  of  algebraically  simple  complex  conjugate 

eigenvalues  iiw^,  u0  f  0,  and  no  other  eigenvalue  is  an  integral  multiple  of  iWQ. 

Moreover,  assume  that  those  eigenvalues  cross  the  imaginary  axis  transversally,  i.e.  if 

X(v  )  denotes  the  branch  of  eigenvalues  for  which  X(V  )  »  iw  , 
o  oc  o 


then 


——  He  X(V  )|„  ,,  f  0.  Then  there  is  an  analytic  branch  of  nontrivial  time-periodic 

av  o  V  *  V 

o  o  oc 

solutions  (u*(£),  v*(E),  p(e),  q(e),  h*<e))  *  s  Y(e),  such  that  Y(e)  e  V  is  a  solution 
of  (2.1)*  -  (2.8)*  for  Vq  •  Vq  (e)  with  temporal  frequency  =  u)  (e ) .  For  e  =  0,  we 
have  VQ  «  voc»  “  ”  “0  and  y  =  0*  This  branch  of  periodic  solutions  is  unique  except 
for  phase  shift  or  changes  of  parametrization. 

If,  at  VQ  =  V^,  all  eigenvalues  other  than  have  negative  real  parts  and  we  have 

Re  X(v  )|v  ^  >0,  then  the  bifurcating  periodic  solutions  are  stable  if 

o  o  oc 

V  (c)  >  V  for  small  c  1  0,  and  unstable  if  V  (e)  <  V 
o  oc  o  oc 

Remark: 

It  is  easy  to  show  higher  regularity  of  the  bifurcating  solutions  by  choosing  function 
spaces  of  higher  regularity  for  the  bifurcation  analysis. 


5.  Reduction  of  the  bifurcation  problem  to  local  form 


In  the  previous  two  sections,  we  have  provided  the  analytical  tools  and  the  estimates 


needed  to  establish  a  bifurcation  theorem.  In  the  following,  we  now  describe  an  algorithm 


for  calculating  approximations  to  this  bifurcating  solution. 

As  before,  we  study  bifurcation  from  plane  Couette  flow,  and  we  consider  the  velocity  of 

the  upper  plate  as  the  bifurcation  parameter.  Plane  Couette  flow  is  the  following 

solution  of  (2.1)  -  (2.8): 

h(x)  -  h  , 
o 


*  n2  Vo  Y 

ui  “  VVVV 


«2  "  0,  p  -  0, 

.  iy  ♦  yyyi  VP 

V1  -  VVVV 


v2  -  0,  q  -  0. 
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For  the  bifurcation  problem,  it  is  convenient  to  introduce  new  variables  representing  the 
perturbation  of  this  solution.  We  therefore  replace  u^  and  v^  by  u1  +  u1 


v  +  v^,  where  u^»  v^  are  given  by  the  formula  (5.1)  in  the  regions 


0  <  y  <  h  (x,t)  and  h(x,t)  <  y  <  1,  respectively.  Moreover,  we  set 


(5.2)  h(x,t)  »  h  +  S(x,t), 
o 


and  6(x,t)  has  zero  mean  value  as  a  function  of  x.  With  this  change  of  variables, 
(2.1)  and  (2.2)  take  the  form 


3u 


3u. 


(5.3)  p  [  u  +  u.  t~  +  e  u,  t —  +  (u*7)  u]  »  n.Au  -  7p,  0  <  y  <  h(x,t), 
1  dx  “x  2  <*y  »  *■ 

div  u  »  0  , 


.  »  3v< 

(5.4)  p  t  v  +  v1  —  +  e^  v2  -J— 1  +  (v*7)  vj  -  n2  Ay  -  Vq,  h(x,t)  <  y  <  1, 


div  v  “  0 

On  the  walls,  we  have 


u  -  0 


(5.5) 


v  »  0 


at  y  -  0  , 
at  y  «  1  . 


In  the  normal  and  shear  stress  conditions,  the  terms  introduced  by  u  ,  v^  are  such  that 


they  cancel.  Moreover,  h'  ■  S',  so  h  can  be  replaced  by  6.  Across  the  interface 
y  -  h(x,t),  we  thus  obtain  the  following  conditions  resulting  from  (2.6)  and  (2.7) 


3u_  3u 


3u_  3u„ 


d  -  A-2)  n,(  aT  +  iT>  *  *»,  5'<  iT  '  ir  > 


(5.6) 


3v  3v  3v  3v 

(1  -  S'2)  72  (  3^  )  *  2n2  > 
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The  spectral  problem  for  the  stability  of  the  null  solution  Is 


p  [Xu+u  j  +  e^  u2u^  ]  +  7p  -  2n1  <Uv  D[u]»0  0  <  y  <  h 

dlv  u  »  0  , 


(6.2)  PtXj+v^  j ~  +  V2V1  1  +  ”  2n2div  D  1^1  “  0  h^  <  y  <  1 

<liv  v  »  0  . 

Here  D  [u_]  Is  the  symmetric  part  of  7u.  On  the  walls  we  have 


(6.3) 

u  * 

0 

at 

y  •  o  , 

• 

V  « 

0 

at 

y  -  t. 

From 

(5.8), 

we  have 

6 

V(Who 

-  kin,  -  v,),  where  k  -  ..  ,  • 

o 


By  inserting  this  into  the  remaining  interface  conditions  and  linearizing,  we  find  the 
following  conditions  at  y  »  hQ: 

(6.4)  u2  -  v2  -  0, 

(6.5)  2?),  D12  tu]  -  2n2  D12(v]  -  0, 


♦  2n jU,j  d12(u*1  -  2n1  u*  d12(u]  +  2n1uJ  D22[u*]}dx 


By  considering  special  forms  of  ji,  v_,  p,  q,  we  find  that  in  0  <  y  <  hQ,  we  have 

—  *  _  »  _ 

(6.9)  pXu*  -  p  u,  t —  +  P  e  u*  u*  -  2n  div  D  [u*]  «  -  7  p», 

—  1  3x  -yll  1  — 

div  u*  “  0 , 
whilst  in  h  <  y  <  1, 

_  .  3v*  _  .  _ 

(6.10)  pXv*  -  p  v1  -g—  +  P  £y  v*  v*  -  2n2  div  D  fv*J  »  -7g*, 

div  v*  ■  0  . 

We  insert  (6.9)  and  (6.10)  back  into  (6.8),  and  compute 

+  /  u  •  ^p*  .  +  /  »*Vq‘  dxdy 

°1  fl2 

L  _ 

”  /  (p»  u2  -  q*  Vj)  dx  . 


This  term  is  added  to  the  right  hand  aide  of  (6.8),  leading  to 
L  _ 

(6.11)  0  -  /  {<-q+2h2  D22  [v])v*  ♦  (p-2h1  D22[u))uJ 

o 

+  (q*  -  2n2  d22(v«])  v2  +  ( -p*  +  2n1  d22  lu*j>u2 

♦  2t)2  v*  d121-1  ‘  ”l  °12[SJ 

♦  2h1  u1  012  [u*]  -  2D2  v1  D12tv*)  }  dx  . 

We  use  (6.6)  to  reduce  the  first  line  of  (6.11),  (6.4)  for  the  second  line  and  (6.5 
the  third  line.  Thus  we  find 
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)  for 


L 

(6.12)  0  -  I  {(-q+2n„  D  (v) )  (v*  -  u*)  -  T  6*u* 

1  2  22  2  2  2 
o 

+  u2<q*  -  2r>2  D^lv*)  -  p*  +  2tl1  D^lu*)) 

+  (v*  -  u*)  2n1  di2  tu] 

*  <Ui  "  v1>  *  2n,  D12t«*l 
+  v1(2n1D12  (u*l  -  2n2  d12[v*])>  dx 

Da  next  writs 

L  L 

J  5"u2*  dx  «  J  6u2*“  dx, 
o  o 

and  Mt  8  -  k(u1  -  v ,  u2  -  (X+uOO  )  k  (u^v^. 

This  leads  to 

L 

0  -  J  «-<I+ 2T12  DJ2  [V]  )(▼•-«•) 

o 

+(u  -v  )  <  -Tk  u  *"  ♦  (X-u(h  )  |-  )  k  (q*  -  2n  D  [v*] 

1  1  2  O  oX  2  22 

-  p*  ♦  2n1  D22fu»J)  ♦  2di  di2  Iu*J  > 

♦  (v*  -  u*)  •  2t)l  012  [u] 

+  vl  ( 2ri ^  D12fi*I  “  2n2  Di2  *  <5x# 

This  yields  the  following  four  conditions  on  y  «  hQ  : 

(6.13)  2n1  d12[u*]  -  2n2  b12[v*1. 
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(6.14)  u*  -  v*. 


(6.15)  u»  -  v». 


(6.16)  -Tk  u*"  +  (X-u(h  )  )  k  (q*  -  2n  D  tv*] 

2  O  <*x  2  22 


-  p  *2 ni  D  [u*)>  +  27^  D12  fu*)  “  0. 


Thus  the  adjoint  problem  is  given  by  the  differential  equations  (6.9),  (6.10),  the 
Dirichlet  conditions  u*  »  0,  v*  »  0  on  the  walls  and  conditions  ( 6. 13) — { 6. 16)  on  the 
interface . 

It  is  easy  to  establish  a  necessary  and  sufficient  condition  for  the  solvability  of  the 
inhomogeneous  problem  corresponding  to  (6.1)-(6.7).  Suppose  that  the  zeros  on  the  right 
of  the  first  equation  in  (6.1)  and  (6.2)  and  on  the  right  of  (6.4)-(6.7)  are  replaced  by 

(6.17)  (x,y),  g2(x,y),  g3  (x),  g^lx) ,  g5(x),  gg(x), 

respectively.  This  inhomogenous  problem  is  solvable  if  and  only  if  the  data  (6.17)  are 
orthogonal  to  the  kernel  of  the  adjoint,  that  is,  when  (6.17)  is  such  that 


(6.18)  J  2  •  u*  dxdy  ♦  J  g  •  v*  dxdy 

°1  fl2 


-  +  J  {g3<5*-2n2  D22t£*)>  ♦  «*  ♦  95  “2 

O 

♦  g  (q#-2n  D  [v#]  -  P*  ♦  2n  D  dx. 

o  2  22  I  22  * 


Me  are  interested  in  the  neighborhood  of  a  critical  point,  where  a  loss  of  stability 


occurs.  There  is  a  critical  plate  velocity  V  -  V  such  that  the  real  part  of 

o  o 


X  vanishes,  and  X  «  iu  .  We  put 

o 


V  -  V  (HR), 
o  o 


so  criticality  is  when  F  «  0.  Ms  get  Hopf  bifurcation  whan  the  loss  of  stability  is 


strict 


„  3X  I  ,  „ 

**  Si  |  r-o  *  °* 

Let 

u  in  ft  (y  <  h  )  , 

C  -  {  "  1 

v  in  fl_  (y  >  h  ) 

2o 

be  a  right  eigenvector  satisfying  (6.1)-(6.7)  at  criticality.  Then  is  the  right 

eigenvector  belong  to  -io)^.  and  are  the  left  eigenvectors  belonging  to 

iio  and  -iui  . 

o  o 


7.  Domain  perturbations  and  Hopf  bifurcation 

He  have  already  demonstrated  that  Couette  flow  can  bifurcate  into  a  time-periodic 
solution  in  which  we  have  travelling  interfacial  waves.  To  compute  this  solution  we 
would,  following  Lindstedt,  sup  the  solution  into  a  fixed  frequency  domain 
( 2w  periodic  in  s) 
udt  ”  ds , 
and  replace 

3u  3v 

(7.1)  (u,  v,  i)  with  u(  > 

in  (5.3),  (5.4)  and  (5.10).  He  then  map  our  problem  into  a  fixed  spatial  doeiain,  using  a 
one  to  one  linear  mapping,  which  takes  boundary  points  into  boundary  points 


(7.2),  y  -  (h(x. 


y  -1  f  hty<  1, 

t>_ho)  hTT  +  yo  <  _  « 

o  ^  ho<yo<1. 


(7.2)2  y  -  (1  + 


h(x,t)-h 


(Ky«h(x,t), 


0<y  <h 
o  o 


He  then  change  variables,  putting  x  ■  xQ  and  y  “  y(xQ,yo) ,  where  y  is  defined  by 
(7.2)  in  equations  ( 5. 3 )— < 5. 10) .  The  form  of  these  equations  changes  under  the  change  of 
variables.  However,  following  ideas  Introduced  by  Joseph  [18],  [19]  we  find  many 
simplifications.  He  shall  now  explain  these  simplifications. 


First  we  introduce  an  amplitude  parameter  which  is  conveniently  set  as  a  projection 


C  “  [  u,  *•]  , 

where 

def  l  2n 

[a,b]  *  2s  /  *  >  ds' 

o 

and  <a^  b  >  are  intergrals  over  both  regions  of  the  type  displayed  in  (6.8).  We  are 
working  in  the  frame  of  $  VIII. 3  of  Iooss-Joseph  [17],  and 
(7.3)  z*  -  eiB 

where  is  the  left  eigenvector  at  criticality  which  was  introduced  at  the  end  of  the 

last  section. 

Hie  bifurcating  solution  may  be  computed  in  a  series  of  powers  of  e .  thus 


(7.4)  m(e) 


(7.5,  vo,e, 


V.  +  .., 


It  follows  from  the  classical  theory  of  Hopf  bifurcation  that  o>  and  VQ  are  even 
function  of  e.  we  have  assusied  this  in  writing  (7.4)  and  (7.5).  Moreover, 


u 

(x,y,s,E) 

uIt)(x  ,y  ,s» 
—  o  o 

V 

(x,y ,s,£) 

•  £t+1 

v[ll(x  ,y  , s ) 

(x,y,s,e) 

(t-M)l 

—  o  o 

m. 

p 

1-0 

p 

o  o 

q 

(x,y ,s,E) 

qCt,(x  ,y  ,s> 

The  functions  of  x  and  y  are  defined  in  deformed  domains  with  a  wavy  interface 
h(x,s,E)  -  hQ  **  6(x,s,e).  The  functions  of  xQ  and  yQ  was  defined  in  the  reference 
domain  with  a  flat  interface  at  yQ  -  h.  Hie  perturbation  of  the  interface  6(x,s,e)  can 
be  eliminated  from  (5.6)  i  that  is, 

4(x,s,e)  -  k  (u1-v1 , (x,s,e) 
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1b  an  identity  for  all  x,s,e. 


The  square  brackets  on  the  left  of  (7.6)  indicate 


differentiation  following  the  mapping  evaluated  at  e  -  0.  For  example. 


fn]  „n 

—  (x  ,y  ,.)  - 

00  de"  e  -  0 


3  u(x  ,y(x  ,y  ,s,c),s,e). 

■  ■  ■  ■  o  00 

3en 

There  is  a  simple  differential  calculus  for  domain  perturbations.  The  partial 
derivatives,  holding  x,  y  fixed,  at  e  -  0  are  also  important.  For  these 


3e 

The  two  types  of  derivatives  are  related  by  the  chain  rule 


^£■0 


[1]  <1>  ^  <1>  - 

u  ”  u  +  y  x-  , 


...  <i>  .2  <o>  ,  <o> 

. .  [2]  <2>  .  -<1>  -  .  .  <1>. 2  *  -  <2>  - 
(7.9)  uJ-u  +  2y  ft  +  (Y  )  *Y  ^ 


u(nl  (x  ,y  ,s)  *  u<n>(x  ,y_,s)  ♦  lower  order  terms, 

—  OO  —  OO 


y 

S<n^x  ,s )  °  ■  ,  h  <  y  <1, 

o  w  o  o 


(7.10)  y<n>(x  ,y  ,s)  -  ( 


4<n^x  ,s)  ,  0  <  y  <  h  . 

on  o  o 

o 


On  the  free  surface  yQ  ■  hQ  of  the  reference  domain ,  we  have 

(7.11)  y<n>(x  ,h  ,s)  -  «<n>(x  ,s). 

o  o  o 

The  equations  governing  the  coefficients  in  (7.6)  are  very  complicated  because  the 
differential  operators  with  derivatives  with  respect  to  x  and  y  in  the  field  equations 


must  be  reexpresaed  by  derivatives  with  respect  to  x_,  y  under  the  change  of  variables 


5} 

I 


/.V.V.aV.V-V/ VN.-V  'A"SVuvji% 
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x,  y  +  xQ,  yQ.  Since  this  mapping  is  invertible  we  can  continue  (7.6)  as 


u  (x,y,s,e) 

(+1 

Utt] ( x,y  (x,y,s,e),s) 

—  O 

v  (x,y,s,e) 

00  c 

J  ((+1)1 

1*0 

[fc] 

v  (x,yQ(x,y,s,e ) ,s) 

etc  . 

etc  . 

Fortunately  it  is  never  necessary  to  solve  the  complicated  equations  which  govern  the 
derivatives  [£]  on  the  right  of  (7.6).  In  fact  we  need  only  to  do  much  simpler 
calculations  for  the  partial  derivatives  <(>.  When  the  partial  derivatives  <£>  are 
known  the  total  derivatives  [t]  may  be  computed  by  the  chain  rule  (7.9).  The  point  of 
simplicity  of  partial  derivatives  is  that  they  don't  perturb  the  operators  which  are 
defined  on  region  and  fi2,  below  and  above  the  free  surface,  see  [18]  and  [19].  For 


examPle'  a  <n>  #  <n> 

3u  3u 

(7.13)  div  u<n>(xo,yo,s)  -  +  — 5— 

o  o 


0, 


whereas 

div  u  lnl (x  ,y  ,s)  +  0  . 

—  O  O 

The  same  type  of  simplication  holds  for  the  perturbation  equations  which  arise  from  (7.1), 
(5.3)  and  (5.4).  For  example. 


3u<2>  .  3u<2> 

(7.14)  P  I-o|=  +  £XU2 


3  <0> 

<2>  *  *  <2>  —  .  <0>  “ ,  <2>  .  <1>„  <1> 

u.  +u.  -  +ex  u2  u'f  +  u  u 

3x 


a  <0> 

-  .  <2>  -  <2> 
+  U2  Is  1  *  T'1  -  7P 


The  unknowns  here  are  u<2>,  P<2>#  u2  and  V2’ 

It  is  not  possible  to  avoid  the  total  derivatives  [t]  in  the  equations  (5.6)— (5.10) 
because  these  are  defined  on  a  manifold,  the  interface,  and  not  in  a  region.  The 
interface  conditions  are  of  the  form 
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(7.15)  g(x,y  =  h(x,s,e),  s,  e)  -  0 

and  the  perturbation  of  y  »  h  with  €  cannot  be  avoided.  In  fact  the  interface 
conditions  are  identities  on  the  interface  so  that  tangential  derivatives  on  them  vanish 
(see  [19]). 

The  following  is  a  recipe  for  perturbations  of  the  domain  in  bifurcation  problems. 
First,  we  introduce  the  frequency  u(e)  into  (5.3),  (5.4)  and  (5.10)  using  (7.1).  We 
then  insert  the  series  (7.4),  (7.5)  and  the  series 


u(x,y,s,e) 
v  (x,y ,s ,£) 
(7.16)  p(x,y,s,e) 

q(x,y,s,e) 

5(x,s,e) 


into  div  u^  -  div  v^  «  0,  (5.3),  (5.4)  and  (5.5)  and  identify  the  perturbation  equations. 

These  equations  hold  in  the  reference  domain.  To  get  the  equations  which  govern  the 

interface  conditions  <  5.6)  —  ( 5. 10)  we  insert  the  series  (7.6)  and  identify.  Then  we 

express  the  derivatives  [i]  with  partial  derivatives  <l>,  using  the  chain  rule.  The 

perturbation  equations  arising  from  interface  conditions  are  thus  defined  on  the  flat 

interface  y  “  h„. 

o 

The  series  on  the  right  of  (7.16)  may  be  expressed  on  the  deformed  domain  by 
inverting  the  mapping,  as  in  (7.12).  In  fact,  the  series  on  the  right  of  (7.16)  is  equal 
to  the  series  on  the  right  of  (7.6),  though  the  partial  sums  of  these  two  series  are  not 
equal  (see  equations  of  [18]). 
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8.  Solvability  of  the  perturbation  equations 


We  must  solve  perturbation  problems  of  the  following  form 

(i)  All  functions  of  s  are  2w  periodic  in  s 

(ii)  All  functions  of  x  =  xQ  are  L  periodic  in  xQ 

(iii)  In  the  region  0  <  y^  <  hQ , 

,  <n>  .  <n> 

3u  3u 

(8-1)  5T-  “°' 

O  O 

3u  <n>  3u<n> 

(8.2)  P[u  t-1  +  u  T~  +  e  u  <n>  u’}  -  n  Au<n>  +  Vp<n> 

o  os  3x  “X  2  I  “ 

o 


=  0  <«  ,  V  ,  x  ,  s). 

I  n  n  o 

(iv)  In  the  region  h  <  y  <  1  we  have  the  same  equations  with 

o  o 

v  (x  ,y  ,s),  v(y  ),  q  and  8  replacing  u  ,  u,  p  and  9,. 
—  o  o  o  z  —  — *1 

.  ,  <n>  -  „  <n> 

(v)  u  «  0  at  y  *  0,  v  »  0  at  y  »  1. 

o  o 

(vi)  On  the  interface  y  «  hQ  +  6 ,  we  have  by  (5.8) 


(8.3)  5  -  k  Iui])'  I°J3*fu1"vr 


We  have  eliminated  5  from  the  interface  equations  (5.6),  (5.7)  and  (5.10)  with 
Of  course 


«<n>  -  k  {[“1<n>]J  °n  V0  “  h0 


(vli)  The  four  interface  conditions  on  yQ  «  hQ  are  of  the  form 

<n>  <n>  0  ,A  _  _ . 

u_  -  v_  »  8  (V  ,  x  ,  s), 

2  2  3  n  o 


T’1D12t-<n>]  -  T12D12[I<n>1  "  W  V  8>' 


<n>  .  _  „  _  <n> ,  <n>  „  _  .  <n>. 

-P  +  n  1 D22  [—  1  +  *  "n2D22(2  1 
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